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Abstract 

We derive the first two moment sum rules of tlie conduction electron retarded self-energy for 
both the Falicov-Kimball model and the Hubbard model coupled to an external spatially uniform 
and time-dependent electric field (this derivation also extends the known nonequilibrium moment 
sum rules for the Green's functions to the third moment). These sum rules are used to further 
test the accuracy of nonequilibrium solutions to the many-body problem; for example, we illustrate 
how well the self-energy sum rules are satisfied for the Falicov-Kimball model in infinite dimensions 
and placed in a uniform electric field turned on at time t = 0. In general, the self-energy sum rules 
are satisfied to a significantly higher accuracy than the Green's functions sum rules. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.45. Gm, 72.20.Ht 



I. INTRODUCTION 



The theoretical description of nonequihbrium strongly correlated electron systems is one 
of the most important problems in condensed matter physics. This problem is not only an 
intellectual challenge, but has the potential for many practical applications. Systems with 
strong electron correlations, like heavy-fermion compounds, manganites, high-temperature 
superconductors and strongly correlated oxide multilayers, demonstrate interesting and un- 
usual properties, some of which have already been applied to electronic and magnetic de- 
vices. Due to the expectation for strong tunability of such systems, they are important 
candidates to be used in modern nanoelectronics, like multilayered structures, quantum 
wires and dots. Some of the properties of these materials can be exploited in spintronic and 
orbitronic devices, where the spin and orbital degrees of freedom are manipulated^. Since 
the size of modern electronic devices can be small (~ 10 — lOOnm), the physical processes 
in these systems can become strongly nonequihbrium because they are exposed to strong 
external fields, which are generated by moderate external potentials (~ 1 V) placed over the 
nanoscale structures. The second consequence of a small system size is that the system will 
have enhanced quantum fluctuations. This makes it difficult to study different properties 
of the system, like transport and optics, since we cannot use phenomenological approaches 
that rely on different relaxation times (Coulomb, phonon, etc.) which are longer than typi- 
cal timescales in the system. Recently, much progress has been made in experimental short 
pulse laser techniques, which allow one to study ultrafast processes in different bulk systems 
and nanostructures. These experiments also need a theoretical interpretation. 

Thus, it is important to have exact nonequihbrium solutions for correlated electron sys- 
tems, which can serve as benchmarks for more general approximation methods. This problem 
is complicated even in the equilibrium case, due to the fact that one needs to treat the kinetic 
energy and the potential Coulomb energy terms in the Hamiltonian on equal footing. The 
simplest models for correlated electrons are the Hubbard model^ and the Falicov-Kimball 
model^ (which is a simplified version of the Hubbard model with localized spin-down elec- 
trons). The equilibrium solutions of these models are known only in the one-dimensional 
case, where an analytical Bethe ansatz approach^ can be used for the Hubbard model and 
in the limit of infinite dimensions, where the dynamical mean-field theory (DMFT) can be 
applied^'^ to both models. 
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Similar to the equilibrium case, much progress in studying nonequilibrium properties of 
correlated electron systems has been made in both cases of low and high dimensions. Differ- 
ent approaches, like perturbation theory, equation of motion and variational wave function 



methods were applied to study the properties of strone 



correlated systems in the case of 



JHIiol). 



Recently, a nonequilibrium 



quantum dot and chain systems (see for example Refs. 
generalization of the Bethe anzatz technique was proposedfii and simulations in one dimen- 
sion with the density matrix renormalization group have been performed^. In the infinite- 
dimensional case, the nonequilibrium properties of the Hubbard^*^ and Falicov-Kimball^ 
models were studied by using second-order perturbation theory in U within DMFT. Recently, 
the Falicov-Kimball model was solved exactly^'i^>i^ii2i>22i2i in the presence of a homogeneous 
time-dependent electric field and in the case of a sudden change in the interaction strength 
U.— In these papers, the nonequilibrium generalization of the DMFT approximation was pro- 
posed, which allows one to obtain the numerical solution of the nonequilibrium problem for 
the Falicov-Kimball model. The numerical method is based on the Kadnoff-Baym-Keldysh 
nonequilibrium Green's function formalism, when the nonequilibrium Green's function is 
defined on the Kadanoff-Baym-Keldysh time contour. We studied different properties of the 
model when a constant electric field is switched on at a particular moment of time. We found 
that Bloch oscillations of the electric current can survive for a long time and develop beats 
with a period depending on the interaction strength; in addition, the Wannier-Stark peaks 
in the density of states can broaden and split, when the Coulomb interaction increases. It 
was also found that the Falicov-Kimball model does not switch from one equilibrium state 
to another when the interaction strength is suddenly changed. 

Since most solutions of strongly correlated problems are numerical, it is important to de- 
velop tests that allow one to check the precision of those solutions. In equilibrium, one of the 
ways to check the accuracy is to calculate the spectral moments of the Green's function^S and 
compare them to exact results. Spectral moments have been used in many different contexts 
than just to test the numerical accuracy of numerical solutions. Harris and Lange^ used 
spectral moments and a projection that forbids double occupancy, to determine properties 
about the spectral moments of the individual Hubbard bands at strong coupling. They also 
determined the equilibrium Green function moments for the Falicov-Kimball model when 
they examined an alloy disorder Hamiltonian. Nolting^ used the spectral moments to de- 
velop different strong-coupling-based approximations to the Green functions of the Hubbard 
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model. This approach has been extended in many different directions to loolc for magnetic 
order or to improve iterated perturbation theory in dynamical mean-field theory when away 
from particle-hole symmetry^i2ii2^i2£i2S. Steven White used the exact expressions for the 
zeroth and the first two spectral moments for the Hubbard model to estimate the accuracy 
of a quantum Monte Carlo solution of the two-dimensional Hubbard model^. Usually, only 
the zeroth and the first two moments have been examined. However, as was argued in 



Refs. 



28 and 



29 



it is also important to know the third spectral moment, since it is con- 
nected with spontaneous magnetic order in correlated systems, and knowledge of the zeroth 
and the first three moments also contain valuable information about the strongly correlated 
bandstructure. The authors of these papers have also established a relation between the 
zeroth and the first moment for the self-energy with the lowest moments for the Green's 
functions. This allowed them to estimate the precision of the solution for the self-energy 
at high energies. Recently, interest in the self-energy spectral moments has been renewed, 
due to an application of these results to the description of experiments on the self-energy of 
high-temperature superconductors arising from angle-resolved photoemission^'^'^. While 
the retarded Green function moments we discuss here are appropriate for the full spectral 
function, the lesser moments (and the greater moments which can be extracted from the 
retarded and lesser moments) are appropriate for photoemission or inverse photoemission 



experiments. The recent work in Ref. 
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examines the lesser moments with a further strong- 
coupling projection that removes doubly occupied states. We do not examine these kinds of 
projections here. Instead we focus on nonequilibrium effects. 

The nonequilibrium case is more complicated than the equilibrium case. In nonequi- 
librium, all Green's functions now depend on two time variables, as opposed to just the 
time difference in equilibrium. Nevertheless, exact expressions have been foundi^ for the 
zeroth and the first two spectral moments of the nonequilibrium lesser and retarded Green's 
functions for Falicov-Kimball and the Hubbard models (coupled to a homogeneous and time- 
dependent electric field). Surprisingly, the retarded moments are time independent for an 
arbitrary time dependence of the electric field. The moments were also used to test the 
accuracy of the nonequilibrium solution to the Falicov-Kimball model in the limit of infinite 
dimensions. However, as mentioned above, it is important to also know the third spectral 
moment, not only to quantitatively improve the measurement of the accuracy of solutions, 
but also to extract information about the quantum state of the system (like the renormalized 
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band structure or the appearance of magnetic order). 



In this contribution, we generahze the results of Ref. |l7l by deriving the third spectral 
moments for the retarded and the lesser Green's functions, and deriving expressions for 
the corresponding zeroth and the first spectral moments of the retarded self-energy for the 
Falicov-Kimball and Hubbard models. Surprisingly, the third-order moment of the retarded 
Green's function (Falicov-Kimball model) and the zeroth (both) and first (Falicov-Kimball 
model) moments of the retarded self-energy remain time-independent. We apply these results 
to benchmark the precision of the DMFT solution of the Falicov-Kimball model in both the 
equilibrium case (at arbitrary doping) and the nonequilibrium case (at half- filling) , when a 
constant electric field is switched on at a particular moment of time. 

The rest of the paper is organized as follows. The equilibrium Falicov-Kimball and Hub- 
bard models and their generalization to include the external electric field are presented in 
Section [III The results for the spectral moments are presented in Sections IIIII (Green's 
functions) andHV] (self-energies). In Section |Vl we give a brief description of the nonequi- 
librium DMFT formalism, present equilibrium and nonequilibrium solutions of the infinite- 
dimensional Falicov-Kimball model and compare results for the moments obtained from the 
numerical solutions with the exact results. Our summary and conclusions are presented in 
Section IVIl 



II. HAMILTONIANS FOR THE MODELS IN EQUILIBRIUM AND IN A UNI- 
FORM FIELD 

The generalized equilibrium Hamiltonian for the spinless Falicov-Kimball and the spin 
one-half Hubbard models can be written in the following unified form: 

^(0) = -Yl t^Mcj - E 4///.- - E ^I^^ - E ///^ + ^ E ///^^l^- (1) 

ij ij i i i 

where in the case of the Hubbard model, the operators Cj (/j) and // (//) correspond to the 
spin-up (spin-down) electron annihilation and creation operators on site i. In this paper, 
we consider the case of a hypercubic lattice, and assume that the electrons can hop to the 
nearest neighbor site. The corresponding hopping matrices are tij = t{j and the chemical 
potentials are /i = /ij for both kinds of electrons (Zeeman splitting can be incorporated by 
choosing different chemical potentials, but for simplicity we keep them equal here). The last 
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term in the Hamiltonian describes the local Coulomb repulsion between spin-up and spin- 
down electrons with a strength equal to U. The Hamiltonian in Eq. ([1]) also corresponds to 
the spinless Falicov-Kimball model, when one sets t{j = 0. In this case, the system consists 
of two kinds of electrons: itinerant c-electrons and localized /-electrons, which locally repel 
each other. In the case of the Falicov-Kimball model, we shall also put /ij = for simplicity, 
since the value of the chemical potential of the localized electrons is not important for the 
spectral moments of c-electrons, which we evaluate below. 

The electric field E(r, t) can be introduced into the Hamiltonian by means of the Peierls 
substitution for the hopping matrices^: 



tij ^ t^j exp 



(2) 



tli exp 



le 



A(r, t)(ir 



(3) 



where the electric vector potential A(r, t) is connected to the electric field in the following 
way: 



E(r,t) = - 



1 aA(r, t) 



(4) 



c dt 

and the scalar potential vanishes. This choice of the electromagnetic potential, when the 
scalar potential is set equal to zero, corresponds to the Hamiltonian gauge. For simplicity, 
we also assume that the electric field is spatially uniform and it lies along the direction of 
the elementary cell diagonal: 



A(r,t) = A(t)(l, !,...,!). 



(5) 



Neglecting the spatial dependence of the vector potential, assumes that we neglect the 
magnetic field effects in the system [since the magnetic field is H(r,t) = V x A(r, t)], 
assuming that the electric field is smooth enough in time, that the transient magnetic field 
can be neglected. This can take place in nanostructures, when an applied external potential 
produces an almost homogeneous electric field due to the small size of the system (see also 



the discussion in Ref. 



The Hamiltonian (in the Schrodinger picture), which describes the electron system cou- 
pled to an external spatially independent electric field, has a rather simple form in the 
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momentum representation (the creation and annihilation operators now create or annihilate 
electrons with definite momentum): 

eA{ty 



k 



k- 

6^ (k 



he 

eA{t) 
he 



P.k,q 



(6) 



where the free electron bandstructures are: 



e k 



eA(t) 
he 



he J 



-2t E^ cos 



a I k 



eA,(t) 
he 



(7) 



d is dimensionality of the system and t is the corresponding hopping parameter. In the case 
of the Falicov-Kimball model, one has to put {k — ^^^j — yU/ = in Eq. ([6]). 



III. SPECTRAL MOMENTS FOR THE GREEN'S FUNCTIONS 

In the case of nonequilibrium, there are two independent Green functions, which describe 
the properties of a many-body system. We use the retarded 

G^{hM) = -i0iti-h) ({ck(ti),4(t2)}) (8) 

and the lesser 

G'<(tl,t2)=^(4(t2)Ck(tl)) (9) 

Green functions as the basis functions. The fermion operators on the right hand side of 
Eqs. ([8]) and iQ are in the Heisenberg representation and the averaging operation (...) is 
performed with respect to the equilibrium Hamiltonian (corresponding to the initial con- 
ditions prior to the field being turned on). It is convenient to use the Green functions in 
Eqs. ([8]) and iQ, since they have important physical interpretations. Namely, the poles of 
the retarded Green function define the energy levels of the system (and thereby determine 
the many-body density of states), and the equal time lesser Green function describes the 
occupation of these levels (and hence determine the distribution function). In equilibrium, 
only one of these functions is independent, since they are connected by a simple relation 
depending on the Fermi-Dirac distribution. 

In order to calculate moments of the spectral functions at different values of time, it is 
convenient to introduce Wigner's time variables for the Green functions in Eqs. ([8]) and ([9]): 
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the average time T = (ti+t2)/2 and the relative time t = ti— 12. The frequency dependence of 
a Green function can be calculated by Fourier transforming the Green function with respect 
to the relative time coordinate, and the time evolution of the function is then described by 
the average time coordinate. In other words, the average time coordinate is associated with 
the physical time in the system. The spectral function for the retarded and the lesser Green 
functions can then be defined in the following way: 

Af'<(T,tu) = ^Im / dte*"*Gf'<(T,t), (10) 

where we have introduced a prefactor t], equal to —1 for the retarded Green function and 
1 for the lesser Green function in order to have positive zeroth moments for both retarded 
and lesser Green functions (see below). The nth spectral moments that correspond to the 
spectral functions in Eq. ( fTOi) are defined to be 

/oo 
dujuj-A^'<iT,u). (11) 
-oo 

It is not difficult to show from Eqs. ( fTOj) and ( JTTl) that there exist the following relations 
that connect the moments with the corresponding Green functions: 

/"OO ;"00 on 

/i^'<(k,r) = ^lmj Jujj Jte^-H-—GZ'<iT,t) (12) 



and 

r „ 1 

(13) 



/zf'<(k,T) = 2^Im 

71 



an 



i=0+ 



(for details, see Ref. llTl ). It is more convenient to use the expression in Eq. (fT2l) for the 
retarded Green function, and in Eq. f|T3|) for the lesser Green function. The time derivatives 
with respect to the operators of the Green functions in Eqs. (j8]) and (Q can be expressed 
by taking commutators of the corresponding fermion operators with the Hamiltonian in the 
Heisenberg picture (the terms proportional to the time derivatives with respect to the theta 
function in the case of the retarded Green function do not contribute to the moments). 
This leads to the following expressions, which connect the zeroth and the first three spectral 
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moments with specific correlation functions: 



/xf(k,T) 
/^f(k,T) 



{{ck{T),cl{T)}) 
1 



(14) 
(15) 



({L^c,(r),4(T)})-({c,(T),L^4(T)}) 
({L^c,(T), cl(T)}) - 2({Lic.(T), Li4(T)}) + ({c.(T), ^^^(T)})] , (16) 
({L=^c,(r), ci(T)}) - 3({L2c,(r), Li4(r)}) + 3({L^c,(r), L2cl(r)}) 



-({c,(r),L=^cI(T)}) 



(17) 



/^o<(k,T) 
/x<(k,r) 

/^2<(k,T) 

/"3(k,T) 



2(4(T)c,(T)), 



(18) 
(19) 

cl{T)L'c,{T)) - 2{(l^{T)) {L\,{T))) + ((l^4(T)) c,(T))l , (20) 



{c\{T)L^c,{T))-{{L'ci{T)]c,{T)) 



{cl{T)L'c,{T))-3{(L'cl{T)] (/.^^(T))) + 3( L^c^T) (L1c,(T))) 



-((^M(^))c.(T)) 



(21) 



wliere = [...[[0, 7Y/f(T)], 7i//(T)]...7-^/^(T)] is tlie multiple commutation operator with 
respect to the Hamiltonian (in the Heisenberg picture), performed n times; the operator 
TiuiT) is given by Eq. ([6]) with all fermionic operators replaced by the Heisenberg-picture 
operators evaluated at time T. The commutation relations can be evaluated directly because 
two fermionic operators at equal times (within the Heisenberg picture) satisfy canonical 
commutation relations. 

Evaluating the commutation and anticommutation operations in Eqs. f lT^ - ([T7|) results in 
the following expressions for the retarded moments: 



/io(k,T) = 1, 

/if(k,T) = [e(k-eA{T))~ii] + Unf, 

/i^(k,T) = [e{k-eA{T))- ^if + 2U[e{k-eA{T))- ^i]nf + U^nf, 



(22) 
(23) 
(24) 
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/if (k,T) = [e{k - eA{T)) - fxf + 3U[e{k - eA{T)) - ^ifnj + 'iU'^[e{k - eA{T)) - 

+ U^Y.[ef{p + q- eA{T)) - 28^ {p - eA{T)) + {p - q - eA{T))] {flU){T) 
p,i 

- U'^Y^ [e^{p + q- eA{T)) -e^{p + q + q'- eA{T)) 

p,<i,<i' 

- e^ip - eA{T)) + e^ip + q' - eA(T))](4V^^^,/,4_,c,+,,)m 

+ U^Yl [^('^ + ^ ~ eyl(T)) - e{k - eA{T)) + 6^{p' - eA{T)) - {p' - g - eA{T)) 
p,p',(i 

+ 2ef{p - eA{T)) - 2ef{p + q - eA(T))] (/J,_,/p'/J+,/p)(T) + U%. (25) 

Summing over momentum yields the following local moments: 

/Xo^(T) = 1, (26) 
/zf(r) = -fx + Uuf, (27) 
^2iT) = Y + /i' - 2[//in; + U^Uf, (28) 

(T) = —{fi - f/n/) + 3Ufmf{fi - f/n/) + 3U^fmf{nf - 1) + f/^ri/ - /i^ 



^^X- (-(///.///.) + (///./]/.) - 2(///.///,) + 2(/;/, ///,)) . (29) 



where i{j = t/exp[— z J^' A(r, t)(ir]. In these equations, we have assumed we are on the 
infinite-dimensional hypercubic lattice, and have evaluated the second moment of the hop- 
ping matrix explicitly; the generalization to finite dimensions is simple to complete (see the 
erratum of Ref . llTl ) . 

As follows from Eqs. (I26 ll -( l28ll . the zeroth and the first two retarded moments remain time 
independent even in the case of an arbitrary external time-dependent field. The third local 
moment [in Eq. (129|) ] is time-independent for the case of the Falicov-Kimball model (t-^ = 0). 
In the case of the Hubbard model, its expression is complex and we cannot immediately tell 
whether they are time dependent (but they most likely are). The last two terms in Eq. (12^ 
are defined by electron correlations and they define the shape of the spectral functions of 
the lower and upper Hubbard bands, the redistribution of the spectral weights between the 
bands and a shift of their centers of gravity^>2E. It is difficult to obtain analytical expressions 
for these terms. 
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In a similar way, one can obtain expressions for the lesser moments from Eqs. (fT8|) - (l2T|) : 



/^o(k,T) 
/i<(k,T) 



2K(T)) 

2[e{k - eA{T)) ~ ^]{nk{T)) 

+UY1 [i4^k+JLfp)iT) + {cUc,fl^J,)iT) 



p,q 



/x<(k,T) = 2[e{k - eA{T)) ~ fi]^n,{T)) 

+ ^U[eik - eA{T)) - /.] $^[(/,V,/p4_,c,)(T) + (/,V,/p4c,+,)] 



(30) 



(31) 



+ \uY,Hk - q - eA{T)) - fi]{fl+jA-,Ck){T) 
p,q 

+ luJ2Hk + q - eA{T)) - ^](/;+,/,4cfc+,)(T) 

-^f/5^[e^(p + g-eA(r)) -e/(p-eA(r))]K4Vp4^+.)(^) " (/pVp4-/fc)m] 
+ ^^' E f(/p+,/p/p+Q/p4-,-QC.)(T) + 2(/,V,/p/^+Q/pcL,c.+Q)(r) 



+ {fp+qfpfp+QfpCkCk+Q+q){T) 



(32) 
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//<(k,r) = 2[e{k - eA{T)) - fA'{clck){T) + 2U{e{k - eA{T)) - f^fJ^i'^kCk+jl+JpKT) 

p,g 

+ 2UY^[s(k - eA(T)) + s(p - eA{T)) - 2// 

p',p 

+ e^ip' - eA{T)) - e^ip' + p- k - eA{Tme{p - eA{T)) - fi){clc,fi,^^_MiT) 

+ 2Uj2Hk - eA(T)) + e{k + q' - eA(T)) - 2ii + (p - eA(T)) 
g',p 

- ef{p + q'- eA{T)W{p - eA{T)) - {p + g' - eA{Tmclc,^,, fl^,U){T) 
+ 2C/2 [e{k - eA{T)) + e{k + q' - eA{T)) - 2^ + - eA{T)) 

- e^ip' + eAiT))]{fi^Jjl,^^J,4c,^,,+,){T) 

+ 2C/2 ^ [£f(^p + q-eA{T))-e^{p + q + q' -eA{T)) 

q',p,q,k' 

- ef{p - eA{T)) + + 5' - e^(T))](4V^^^,/,4c,+,,4,_^c,0(T) 

+ 2C/2 + q + q'- eA{T)) - - {p' + q' - eA{T)) + {p' - eA{T)) 

p',g',p,g 

_ e^ip + q - eA{T)) + ef{p - e>l(r))](4V^,/,,/,V/,4c,+,+,,)(r) 

+ 2U' Y ^fp+QfpfW,'fp'fUM<^Q+k+g+g')^T). (33) 
p' ,q.' ,p,(i,P,Q 

The corresponding local lesser moments are 

/i<(T) = 2ne(T), (34) 
/x<(r) = -2j]t,,(cfe)-2/.ne + 2C/^(///4c,), (35) 

t,<{T) = 2^t,;t,,(clc,) +4/. J]f,,(clc,) + 2/xV(T) -2t/ J] \^tij{fl f4c^) +i,,{fl Uc]c,) 

+2U{U-2i^)Y{fIfi4ci): (36) 
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/x<(T)= - 2 J2 iatiJ^Ach){T)-QliY.tdiMc,){T)-Q^^ 

i,j,l,m i,j,l i,j i 



+ 6UfiJ2u,{clc,fjf,){T) + 6f//xE t.,(clc,///.)(T) - 2U'Y.^,,{fjf4c,){T) 



- 2f/2 5^t,,(///,clc,)(T) + 2UY,UM4cif]fj){T) + 2f/5^t,,t,,(clc,///.)(T) 

- 2U^Y.~^,,{f]f,flf4c,){T) 

+ 2Uj2uMi{4c,fhi)(T) - 2UY,Wiu{c\c,f}h){T) + 2U ^^Y,i{M^^fl f))iT) 

- 2U^ J2 (4^/]/.) (T) + 2t/ tl^X- (^1^^///^) (^) + E (S^c.///.) (^) 



where iij = tijexp[—i j^' A(r,t)(ir]. 

Contrary to the case of the retarded moments, even the zeroth and the first two local 
lesser moments in Eqs. flMl) - fl5Bl) cannot be expressed solely in terms of the model param- 
eters, and they depend on different correlation functions. Therefore, in order to check the 
accuracy of calculations in the lesser case, one can only compare the numerical results for 
the lesser moments obtained by direct calculations by using Eq. ffTTl) with the correspond- 
ing numerical results obtained by the evaluation of the Green function time derivatives in 
Eq. f[T^ . However, the results in Eqs. flMl) - fl?n) still contain practical importance because 
they provide a simple way to calculate combinations of different correlation functions. The 
reason for this is due to the fact that the correlation functions on the right hand side of 
Eqs. f l34l) - fl37j) can be expressed in terms of the local lesser Green functions and their time 
derivatives by using Eq. (fT2|) . the equation of motion and/or the Dyson equations for the 
Green functions [see Eqs. fl5^ . fHl|) - fH^ below]. For example, as shown in Ref . we can 
connect the average potential energy with the Green functions and self-energies: 



. d /, eA(ii; 



t2=tl 



(38) 



(37) 
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a generalization of the well-known equilibrium result. 



IV. SPECTRAL MOMENTS FOR THE RETARDED SELF-ENERGY 

It is possible to derive expressions for the lowest retarded self-energy moments, by using 
the Dyson equation, which connects the retarded Green function and self-energy, and the 
results for the retarded Green function moments derived in the previous Section. 

In order to derive the nonequilibrium Dyson equation for the retarded Green function, it is 
convenient to write down the Dyson equation for the contour-ordered lattice Green function 
in the Larkin-Ovchinnikov representation, where all the time arguments are defined on the 
real branch of the time contour: 

/oo i*oo 
dh / dhGl{h,h)%{h,U)G],{h,t2) (39) 
-oo J — oo 

and all the Green functions and self-energy functions are 2x2 matrices 

G.{t,M)=\ ^-'(^-^^) I, (40) 

Gt{UM) 

Sk(ti,t2)= ' / h (41 

with matrix elements which consist of the retarded, advanced 

Gi^(ti,t2) =^^(t2-ti)({ckH(ti),cU(t2)}), (42) 

and the Keldysh 

h) = -z ( [c^M, 4h(^2)] ) (43) 
components (and similarly for the self-energy). The function G^'^ in Eq. fl5^ is the electron 
Green function in the noninteracting case {U = 0, but with 7^ for the nonequilibrium 
case). The expression for this function can be obtained analytically (see, for example, 
Refs. y and Q). 

The nonzero matrix components of the Dyson equation (139|) can be written in the follow- 
ing form: 

^^(^1,^2) = (ti,t2) + [Gf EM(ii,i2), (44) 
G^{h,t,) = Gf{h,t,) + [GfE^G^]{t^,t,), (45) 

Gi^ihM = [1 + [1 + ^^G^Ktuh) + [G^S^Gi^](ti,t2) (46) 
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where we suppressed integrations over internal time variables implied by the continuous 
matrix operator multiplications. 

In order to find the retarded self-energy spectral moments, one only needs Eq. fl44p . It is 
convenient to rewrite this equation in a combined frequency-average time representation 

G^iT,uj) = G^%T,uj) + j dT j di j dn j due-'^'e'"^ 



(47) 



where we restored the internal time/frequency integrations. 

Similar to the equilibrium cas^^^^^, one can expand the Green functions and the self- 
energies at large values of the frequency u in terms of the corresponding moments: 

Gfp-.-) = t^^^' («) 

m=0 

(T, u) = S^(T, ^ = 00) + ^ (49) 

where the moments /i^(k, T) and C^(k, T) correspond to the retarded Green function and 
self-energy in Eq. (1441) . In particular, we have 

C^(k,T) = Im / (ite^^*cj"S^(T,t). (50) 

J-00 

The large-a; expansions in Eqs. fHHl) and fj49l) can be obtained by using the following spectral 
identities (valid for retarded functions that decay rapidly enough for large relative time): 

G«(r..)^-ir..'l^^^^ai^, (51) 



S^(T, ..) = -i r duj' ^'^^^^^f^ + E^(r, = 00), (52) 

J-00 to -Uj' 

where we take u large enough that the Green's function and self-energy on the 1. h. s. are 
real. In fact, by making expansions in powers of (l/cu) on the right hand sides of Eqs. (15T|) 
and (l52l) and by using the moment definitions in Eqs. ( fTTj) and ( fT2l) . one can obtain the 
expansions in Eqs. (HHll and ( H9ll . The self-energy expansion in Eq. (H9ll contains a frequency- 
independent term S^ (T, u = 00), which corresponds to the mean-field term of the self-energy 
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[see Eq. ( 158|) below]; this form arises because the self-energy generically approaches a real 
constant nonzero value as \uj\ —>■ oo. 

Then, one can insert these expansions into Eq. fl47p and consider separately the terms, 
which have the same order in (l/cu). In order to do this, it is necessary to expand all the 
functions under the integrals in powers of (l/cu). For example. 



m=0 ^ 2/ 



(1 + (fi - D/o;; 



m+1 



m=0 

m=0 ^ ^ 

To calculate the frequency-independent term and the zeroth and the first spectral mo- 
ments for the retarded self-energy, it is necessary to make an expansion of the functions 
in powers of l/u in Eq. fj471) up to fourth order. All the time and frequency integrals in 
Eq. (147|) can be easily performed, and we get the following equations which connect the 
Green functions and self-energy spectral moments: 

/.«(k,T) = /i^(k,T), (54) 
/xf (k, T) = /if (k, T) + /i^(k, r)S^(T, CO = oo)/xo^(k, T), (55) 
/xf (k, T) = /If (k, T) + /if (k, T)Sf (T, = oo)/if (k, T) + /if (k, T)Co^(k, r)/if (k, T) 

+/if (k, T)Sf (T, a; = oo)/if (k, T), (56) 
/if (k, T) = /if (k, T) + /if (k, T)Sf (T, ^ = oo)/if (k, T) + /if (k, T)Co«(k, T)/.f (k, T) 
+/if (k, T)Cf (k, T)/if (k, T) + /if (k, T)Sf (T, a; = oo)/if (k, T) 
+/if (k, r)Co^(k, r)/if (k, T) + /if (k, r)Sf (T, u; = oo)/if (k, T), (57) 

where the matrix /if(k, T) is the nth spectral moment of the retarded Green function in 
the noninteracting case. One can straightforwardly derive expressions for the retarded self- 
energy moments from Eqs. ( l5^ -( l57|) by using the results in Eqs. (I24l)-(l25l) for the retarded 
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Green function moments. After some long algebra, we find 

S^(T,tu = oo) = Uuf, (58) 

Co^(k,T) = nf{l-nf)U\ (59) 

Cf(k,r) = U^nf{l-nf)[U{l-nf)-^^] 

+ 5^[e^(j9 + q - eA{T)) - 2e^(p - eA{T)) + {p - q - eA{T))]{fif,){T) 
p,q 

- U^Yl [^^(^ + ^ ~ ^^^^)) - £^(p + g + g' - eA(T)) 

- £^(p - eA(r)) + e^ip + g' - eA(T))](4V^^^,/p4_^Cfc+,,)(r) 

+ ^ [e{k + q- eA{T)) - e{k - eA{T)) + e^ip' - eA{T)) - {p' -q- eA{T)) 

+ 2ef{p-eAiT))-2e^{p + q-eA{T))]{jyj,>fl^J,){T). (60) 

The expressions for the local moments are 

C^{T) = nf{l-nf)U'. (61) 
Cf(T) = U'nf{l-nf)[Uil-nf)-f^] 

ij ij 



It is worthwhile to notice that the local retarded self-energy moments are time independent 
(except for the first moment in the case of the Hubbard model, for which we are not sure 
about the time dependence). This may be a surprising result for the Hubbard model, since 
the second-order perturbation theory is frequency- dependent, but the total weight of the 
self-energy reamins constant and depends just on the electron densities and the interaction. 
Other interesting observations are that the mean-field term S^(T, = oo) is equal to the 
first order (Hartree-Fock) term of the self-energy in the expansion in U, and that the ze- 
roth moment corresponds to the zeroth moment of the imaginary part of S^(co',T) in the 
truncated second-order perturbation expansion^ (for the Falicov-Kimball model case). This 



is in agreement with a result of Ref. 



36 



% where it was shown that in equilibrium the exact 
coefficient of the term proportional to 1/uJn in the large Matsubara frequency expansion 
of the electron self-energy of the Hubbard model can be obtained from the second-order 



skeleton diagram for the exact Green function. Finally, it was shown in Ref. 
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37, in the 



insulating phase, that the imaginary part of the d ^ oo equihbrium retarded self-energy 
acquires an additional term proportional to 6{uj) in the frequency representation (at half 
filling; away from half filling a delta function appears but not at = 0). In particular, 
in the case of the Falicov-Kimball model, the weight of the delta function term is equal 
to —7i[U^7if{l — Tif) — 1/2] and it produces a term that requires special care to include in 
the zeroth self-energy moment, when one performs the integration over frequency of the 
self-energy. Note that the delta function implies that the finite-frequency integration of the 
zeroth self-energy moment remains fixed at 0.5 for the Falicov-Kimball model in the insu- 
lating phase at half-filling, and all of the additional spectral weight comes from the delta 
function piece at = 0. Away from half filling the delta function typically contributes to all 
moments because it appears at a finite frequency. In the nonequilibrium case, the situation 
is more complicated, because we cannot prove that such a term is also present in this case. 
To see whether such a term is present, one needs to examine the large relative-time limit 
of the nonequilibrium retarded self-energy, which would have a constant term equal to the 
weight of the delta function when the delta function appears at = (and would be a term 
proportional to exp iujt when the delta function is at a finite frequency), but we do find good 
overall agreement for the sum rules, so this issue is not important in verifying the accuracy 
(when one performs calculations in the time representation). 

Unfortunately, it is impossible to derive analogous expressions for the lesser self-energy 
spectral moments 



since in this case the expansions similar to Eqs. ( l48l) and ( l49l) do not exist. In fact, the 
representations in Eqs. ( l5Ti) and ( l52i) . which lead to Eqs. f HHl) and ( l49l) . are not valid in 
the cases of the lesser Green function and self-energy, because the lesser functions are pure 
imaginary and hence not analytic. Note that we could try to define an auxiliary Green's 
function that has the imaginary part of the lesser Green's function and a real part determined 
by the spectral function defined by the integral of the imaginary part, but doing so does not 
produce any new results for the spectral moments of the lesser self-energy. 

An alternate approach is to express the lesser self-energy in terms of the retarded Green 
function and self-energy by using the system of Dyson equations in Eqs. ( H4|) -(H6l) and the 
equation which connects the lesser Green function with the retarded, advanced and Keldysh 




(63) 
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Green functions, 



(64) 



and then try to express the lesser self-energy moments in terms of moments for the retarded 
and lesser Green function and the retarded self-energy. In this case, one can find the following 
Dyson equation for the lesser self-energy: 



G<ih,t2) = [1 + [1 + J:^G^]ih,h) + [G^i:<G^]ih,h 



^Rsr</^A] 



Using the equations of motion for the Green functions, 



d 



S{t - h) ( I— - e{k - Aiti)) + /i ) - tij 



G^{ti,t2) 



z J- + e{k - A{t2)) - /i ) 6{t2 - t') + Si^(t2, t') 



- 1') 



d 



S{t - h) { I— - e{k - A{t,)) - ) - S<(t, ti, 



(65) 

(66) 
(67) 
(68) 



one can get the following formal expression for the lesser self-energy: 



^^ih,t2) 



X 



S{h - t) [i^^ - e{k - A{t)) + ) - S^(ti, t) 
+ e{k - A{t')) - fi ) 6it' - t2) + Si^(t', t2 



G<it,t') 



(69) 



Using this result, one can calculate the lesser self-energy moments similar to what was done 
for the Green functions: 



C<(A;,T) = -Im 

TV 



■s^(r,t) 



(70) 



Jt=o+ 



where T and t are the average and the relative time coordinates. 

Unfortunately, this approach also does not provide any useful results for the self-energy 
moments. In fact, even in the equilibrium case, one finds from Eqs. (169|) and (1701) the 
following trivial result: 



(71) 



[In order to obtain this expression, one needs to use the following equilibrium relations: 
Ej^(a;) = S^*(w) and G<{uo) = if{uj)lmG^{uj)]. The result in Eq. ^ can also be obtained 
directly from the equilibrium relation Tj^{uj) = if{uj)lm'E^{uj). Unfortunately, it is impos- 
sible to get analytical results for the lesser self-energy moments from Eq. (I7T1) . except in 
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the high-temperature hmit, when they can be expressed in terms of the retarded self-energy 
moments [via a series expansion for f{uj)]. 

Since the exact analytical results for the lesser moments cannot be found even in the 
equilibrium case, one can try to make some approximations in order to obtain them. The 
standard approximation for the lesser Green function is the generalized Kadanoff-Baym 
(GKB) approximation^^: 

^^(^1,^2) = -^ [G''{h,h)G<{h,t2) - G<{h,h)Gi{h,h)] . (72) 



Substitution of this result into Eq. (!69|) and using the equations of motion in Eqs. (I66l) - (l68|) 
gives the following approximate result for the lesser self-energy: 

E<(ti,t2) = -^ [S«(ti,t2)G<(t2,t2) -G<(ti,ti)Si^(ti,t2)] -2z5(ti (73) 

or in the frequency-average time representation: 

S<(T,c.) = 2amS^(T,c.)nk(r) - 2z^^. (74) 

After summation over momentum the last term disappears, due to conservation of the total 
particle number, therefore in this case 

C<{T)=2j2c^{k,T)n^{T). (75) 

k 

Since the zeroth and the first retarded self-energy moments are momentum-independent, 
one can obtain the following GKB result for the corresponding lesser moments 

C<(T) = 2C«(T)n,. (76) 

The GKB approximation gives good results for the Green's function moments in the case of 
weakly interacting systems. Therefore, the relation Eq. fl76|) should be approximately valid 
in this case. There is one subtle issue with regards to the GKB and DMFT. In DMFT, 
the self-energy is local, and hence momentum independent. But the GKB approximation 
to the self-energy in Eq. ([71]) appears to be momentum dependent. Hence, it is not clear 
how accurate the local self-energy moments will be within this approximation, but because 
the GKB corresponds to a mean-field-like decoupling of correlation functions for the Green 
function moments-^, it is possible that the approximation remains reasonable for the local 
self-energy, at least for weak couphng. 
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Thus, generally speaking, similar to lesser Green function moment case, one cannot obtain 
analytical expressions for the lesser self-energy moments. Moreover, it is even impossible 
to express these moments in terms of correlation functions. Hence, in order to check the 
accuracy of the numerical calculations, one can only compare the numerical results for the 
moments with the numerical evaluation of the self-energy time derivatives in Eq. (1701) . which 
is not a stringent test. 



V. NUMERICAL RESULTS FOR THE FALICOV-KIMBALL MODEL IN INFI- 
NITE DIMENSIONS 

In this Section, we shall use results for the local moments obtained in Sections IlimiVI 
to check the accuracy of the equilibrium and nonequilibrium numerical solutions of the 
Falicov-Kimball model in the limit of infinite dimensions. In this limit, the electron self- 
energy is locaF^, which allows one to solve the problem numerically in both equilibriunt^ 
and nonequilibrium case9^ii^>i^ii^i22i2i. The case of infinite dimensions is important, since 
many physical properties of the model are qualitatively similar as in the 2D and 3D cases 
(see, for example, Ref. I^. 

"^max ^max 



"'max"'P 

FIG. 1: The complex Kadanoff-Baym-Keldysh contour for the two-thne Green's functions in the 
nonequilibrium case. 

In order to study the time-dependent properties of the model in infinite dimensions, 
one needs to solve a generalized system of nonequilibrium DMFT equations for the contour 
ordered Green's function G{ti, ^2), self-energy S(ti, ^2) and an effective dynamical mean-field 
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A(ti,t2): 



k 

^0(^1,^2) = [G-' + ^-\tut2), 

^(^1,^2) = {l-wi)Go{h,t2;fi)+wi[Gglp{fi-U)-\]-\h,h), 



(77) 

(78) 
(79) 
(80) 



where all time arguments are defined on the complex Kadanoff-Baym-Keldysh time contour 
(see Fig. [T]). On this contour, the time increases from the top left point (— t^ax) along 
the contour to the bottom point of the imaginary axis {—tmax — iP)- In Eqs. fl77l) - (|5Ull . 
G^\ti^t2) is the noninteracting electron Green's function in the presence of an external 



field and Goimpiti,t2] fi) is the free impurity Green function; is a chemical 
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potential and 
andbll). 



Wi is the average number of the /-electrons per site (for details, see Refs. 

As mentioned in Section [Tll we shall consider the case of a spatially uniform electric field 
directed along the elementary cell diagonal, as in Eq. We also assume that the system 
starts in equilibrium with an inverse temperature (3 and then a constant electric field is 
turned on at time t = 0. 

In the case of an external field, as given in Eq. ([5]), the free electron spectrum [in Eq. ([7j)] 
has a simple momentum dependence: 

eA{t) 



k- 



hc 



cos 



ea 



sm 



ea 



where 



and 



e(k) = -2t^cos(aA;' 



e(k) = -2t^sin(aA;' 



(81) 



(82) 



(83) 



are two energy functions. It is possible to show that in the case of an infinite dimensional 
hypercubic lattice, the joint density of states for these two energy functions has the following 
formi^: 

1 r ^2 ^21 

(84) 



exp 



e 

^2 



t 



*2 



where t* is a scaled hopping parameter, connected with the hopping t in the Hamiltonian 
Eq. ([1]) as t = t* /2\/d. The momentum summation in Eq. (I77|l can be performed by using 
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the joint density of states Eq. (18^ : ^^-^k = J de J dep2{e,e)F^ g, since in our case the 
noninteracting Green's function on the r h s. of Eq. fl77|) has simple momentum-dependence, 
which can be expressed in terms of the two energy functions in Eqs. fl82p and (1831) . The 
energy integration can be performed by using Gaussian integration^^'^^. We typically use 
about 100 points per dimension. 

In addition, one needs to choose the proper discretization of the time contour Fig. [1] The 
results depend strongly on the discretization step when the step size is not small enough. 
Choosing a given discretization and a tmax determines the matrix size for the given calcula- 
tions. We typically work with general complex matrixes of size 900 x 900 up to 5700 x 5700. 

A. Equilibrium case 

First, we consider the equilibrium case, when there is no external field. In this case, the 
system of equations ( 177l) -( l80l) reduces to the equilibrium DMFT equations^ with no average 
time dependence, so functions of two time arguments can be replaced by corresponding 
functions of one frequency, F{ti,t2) F{uj). The numerics are under good control and 
one can obtain quite accurate solutions. The most important numerical checks that can be 
performed arise from a comparison of the spectral moments calculated directly by integrating 
the real-frequency solutions, with results for the moments that can be determined exactly 
via parameters of the model for the retarded moments or by an evaluation of the relevant 
correlation functions using a Matsubara frequency formalism for the lesser moments. 

Now we show how to calculate the required correlation functions in Eqs. flMl) - fl37j) using 
the Matsubara Green's functions. One starts from the imaginary time-ordered Green's 
functions 

G,,{r) = -{TMr)4m, (85) 

where the imaginary time-dependent operators satisfy Cj(r) = e^'^Ci{0)e~^'^ according to the 
Heisenberg representation. Because these functions are antiperiodic on the interval [0,/3], 
we employ a Fourier expansion in terms of the Matsubara frequencies: 

Gijir) =Tj2e''""" [ dke-'^^^'-'^^^Gki^uJn), (86) 

where ujn = 7rT(2n+l) is the fermion Matsubara frequency. Here, the momentum-dependent 
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Matsubara Green's function satisfies 



1 



(87) 



zcj„ + /i - e(k) - Sk(icj„) ' 

and in DMFT the self-energy has no momentum dependence. 

We start by deriving the equation of motion for the Green function in Eq. fl85|) and 
extracting the expression for the local four-operator correlation function by evaluating the 
Green function at r = and removing the single-particle terms: 



T 



(88) 



n.k 



The correlation functions for operators on different sites, like {fi ficlcj), can be found by 
introducing an extra term — hiflfi with a local field hi into the equilibrium Hamiltonian 
and then evaluating derivatives with respect to hi and taking the limit hi 0. For example, 
straightforward algebra shows that 



0' 



(89) 



where (w,) = (///i) = n/, (see Refs. 
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41 



and 



421 and the Appendix for details). Using 



these identities allows us to find explicit expressions for all of the relevant correlation func- 
tions using Green functions and self-energies determined at the Matsubara frequencies. We 
present the final results for the case of the Falicov-Kimball model in infinite dimensions. 
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where the self-energy is momentum-independent: 



ric = T^G'k(2 

n,k 

^tij{4cj) = -T^£(k)Gk(ia;„), 

i,j ri,k 



71, k 



i,l,m,j n,k 
i 



n,k 



1,3 



T 



-jj$^S(^a;„)£(k)G'k( 



(90) 
(91) 
(92) 
(93) 

(94) 
(95) 



n,k 



^Ujtji{c\cif]f^) 



T 
U 



n,k 



T 



(96) 
(97) 



n,k 



We next perform the momentum summation in Eqs. (I90l) - (l97l) to express the results in terms 
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of local quantities: 



ric = T^Gn, 

n 

^ tij{clcj) = T ^ [1 - {iUn + 11- Sn)Gn] , 
i,j n 

tiltij{c\cj) = -T^{iUn + /i - S„) [1 - {iUn + fi- , 



(98) 
(99) 
(100) 

(101) 

(102) 



T 



- E [1 - (^^n + - S^)^^] , (103) 



T 



= E ^"(^'^^ + // - S„) [1 - (za;„ + // - E„)G„] , 



T 
IP 



(104) 
(105) 



where G„ = ^k^k(^'^n) and E„ = E(ia;„). These expressions can then be employed to 
efficiently determine the lesser moments from an independent Matsubara frequency calcula- 
tion. 

We find, for all cases that we consider, all of the different Green's function and self- 
energy moment sum rules are satisfied to essentially as high an accuracy as we want (the 
delta function contributions to the self-energy moments must be included to get the correct 
answer; this becomes complicated for particle-hole asymmetric cases when U is large enough 
for the self-energy to have developed a pole because one needs to accurately determine the 
location and weight of the pole to obtain the correct sum rules). In some cases, we need 
to use many Matsubara frequencies in the summations to achieve sufficient accuracy, or we 
need to have a small frequency grid spacing for the real-frequency Green's functions. The 
sum rules hold in the case of half-filling and away from particle-hole symmetry and they 
hold equally well for metallic and insulating cases. 
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B. Nonequilibrium case 



In this Subsection, we compare the numerical results for the moments (at half filling) 
with exact analytical results obtained in the case when a constant electric field is turned on 
at time t = 0. Since we calculate the contour-ordered self-energy, we need to extract the 
correct retarded quantities to compare with the moments that do not depend on correlation 
functions (which we have no independent way to evaluate). This is simple to do for the 
Green's functions. For the self-energies care is needed. The constant term in the self-energy 
in the frequency representation becomes an equal time delta function in the time formalism. 
The zeroth moment corresponds to the equal-time retarded self-energy (most easily found 
by taking the difference of the greater and lesser self-energies) and the first moment is found 
from the first derivative. We need to evaluate the derivative carefully, because we need to 
remove the delta-function piece first. We handle this instead by using linear extrapolation 
from finite relative times to the vanishing relative-time limit, so we do not need the data at 
equal times to find the derivative. More sophisticated techniques would be needed to find 
the higher moments, but we don't need those here. 

In general, the self-energy moments are satisfied to very high accuracy, even if the step 
size is large. Errors are often less than 0.1%, which is much lower than what one finds for 
the Green's function moments (where we often need to work hard to get errors below the 
1% leveU^ii^i^) . We can extrapolate the results to the limit At 0, which produces even 
higher accuracy. The results are most accurate for the constant piece to the self-energy. 
Then the zeroth moment, and finally the first moment. But the results of our investigations 
indicate that the Green's function moments are a much more accurate test of the accuracy 
of the solutions than the self-energy moments. While we could show similar scaling plots as 
were created for the Green's function moments^^, it does not seem to be necessary because 
the improved accuracy is so much better for the self-energies that one does not learn too 
much from such an exercise. 

VI. CONCLUSIONS 

In this work we have shown how to extend the Green's function moment sum rules to 
third order for both the Hubbard and the Falicov-Kimball models and used these moments 
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to examine the retarded self-energy moments through first order. Our analysis holds both 
for equilibrium and nonequilibrium situations. We find for the Falicov-Kimball model that 
the moment sum rules remain time independent in nonequilibrium, which is a surprising 
result. In the case of the Hubbard model, it appears that the third order moments will 
be time dependent, but we cannot explicitly confirm this. When we compare the sum 
rules to numerical calculations for the Falicov-Kimball model with DMFT, we find excellent 
agreement both in equilibrium and in nonequilibrium. In fact, the Green's function sum 
rules are a much better indicator of overall accuracy than the self-energy sum rules. 

The sum rules are only relevant for quantitative comparisons of retarded functions. In the 
case of lesser functions, we are able to make comparisons of the Green's function sum rules 
to the relevant correlation functions evaluated with a Matsubara frequency formalism when 
the system is in equilibrium, but we cannot extend that approach to the nonequilibrium case. 
We are unable, even in equilibrium, to find any useful sum rules for the lesser self-energy. 
Instead we find just trivial relationships that arise from the definitions of these quantities 
(which are well known in equilibrium and unknown in nonequilibrium). 

In the future, we will examine how these sum rules can be extended to inhomogeneous 
situations, with relevance to inhomogeneous DMFT (and other techniques) as applied to 
mutlilayered nanostructures or ultracold atomic systems in a harmonic trap. In addition, 
utilizing these sum rules can allow one to obtain more accurate results for the high-frequency 
hmit of the Green's functions, self-energies and dynamical mean fields. We will illustrate 
this use in another publication, which allows one to employ a minimal number of Matsubara 
frequencies yet maintain high accuracy of solutions. 
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APPENDIX A: CALCULATION OF EQUILIBRIUM CORRELATION FUNC- 
TIONS USING THE MATSUBARA FREQUENCY FORMALISM 



In this Appendix, we present details of the derivation of the correlation functions in 
Eqs. flM|) - flW]) . The expression in Eq. fl55]) for the first correlation function in Eq. flM|) can 
also be determined by introducing a fictitious field — ^ ■ hif] into the Hamiltonian and 
taking derivatives with respect to hi and then setting all hi = 0. Each derivative with 
respect to an hi brings down an operator fjf- into the operator average (plus a correction 
term when the derivative acts on the partition function). This approach is more general 
than the equation of motion approach used to derive Eq. fISS]) and will allow us to derive 
expressions for the other correlation functions in Eqs. (pH]) -f l97p . As shown in Ref. |42|, a 
correlation function that contains a product of two c-electron operators and one /-electron 
number operator can be expressed in terms of a derivative of the c-electron Green function 
with respect to the fictitious field [for example, Eq. (l89l) ]. 

In order to explicitly calculate the fictitious field derivative of the Green function, one uses 
the standard trick of writing G = GG^^G so that derivatives of G are replaced by derivatives 



of G~^ which involves a derivative of the self-energy (see Ref. |42| ). Because we have added 
the fictitious fields to the Hamiltonian, and they are not translationally invariant, we lose 
translational invariance in the system prior to taking the derivatives (it is restored once we 
set /ij = 0). Hence, we need to work in real space rather than momentum space, and we 
need to allow the dynamical mean fields and the self-energies to have a site dependence. 
This implies that we can write the local Green function at site i via 

Gii{lUJn) = -. — , ,^ . . (Al) 

in the Matsubara frequency representation. 

Now consider the case where we add an /i-field only at site i. Since the h field will modify 
Hfi, the Green function and self-energy at site i are changed by hi. What about the Green 
function and self-energy on neighboring sites? Using the Dyson equation, one can show that 
the change in the Green function at site j, 6Gjj{iu!n), is equal to 

6Gjj{iuJn) = Gji{iuJn)\h,=QST.i{iuJn)Gij{iuJn)\h,=0- (A2) 

But Gij is proportional to the hopping t raised to the power equal to smallest number of 
hops between site i and site j. So, for example, if j is a nearest-neighbor of site i, the 
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right hand side of Eq. (lA2p is proportional to = t*'^ /Ad — > as c? ^ oo. Hence, we learn 
that 5Gjj{iuJn) = for j 7^ i and large dimensions. If Gjj is unchanged, then is also 
unchanged. This means that 

(A3) 



d'Lj{iuJn) . 



We now show how to derive one of the off-diagonal c-f correlation functions. We want to 
calculate 



I3dh, 



+ nfi 



Gji{iuJn), 



(A4) 



tj n ij 

which follows directly from the definition of the operator average and an explicit computation 
of the derivative (the term multiplied by ri/j arises from the derivative of the partition 
function). Now we focus on the derivative term, and use the GG^^G trick 



d 



I3dh, 



■GjiiiUn) 



d 



I3dh, 



Gji{iuj^ 



Gjk{iuJn)Gf^,^ {iuJn)Gii{iuj„ 
d 



kl 



I3dk 



Gii{iuJn), 



(A5) 
(A6) 



where we used the fact that the derivative of the self-energy was nonzero only for k = I = i. 
Since the self-energy is an implicit function of Ga and n^j one can compute the derivative of 
the self-energy with respect to the field by using the chain rule and re-expressing in terms of 
derivatives of the self-energy with respect to the Green function and the /-electron filling. 
The algebra is quite long and is contained in Ref. 

d 
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Gii{iuJn) 



/3dk, 



The end result is that 



u 



(A7) 



Plugging this result into Eq. (]A6p . and then converting the summation over i and j to a 
summation over momentum, produces Eq. ( I95|) . 

The only equation that requires some more formal development is Eq. fl97|) because it 
involves two /-electron density operators, and hence derivatives with respect to two h fields. 
Using the fictitious fields, one can immediately show that 



d 



pdh. 



d 



0" 



(A8) 



All the terms in this expression, except the term proportional d'^Gji/ dhjdhi, can be expressed 
in terms of the Green function and self-energies by using the results above. In order to find 
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the second derivative of the Green function, one can show (similar to the case of the first 
derivative) that in the hmit of infinite dimensions: 



d'^Gji{iuj. 



dhjdhi 



\h=0 



_d_ 

dhi 



dhi 



\h=0 



dhj dhi 



dhjdhj 



\h=0- 



(A9) 



The last term in this equation is equal to zero, since the second derivative of the self-energy 
d'^T,i{iuJn) / dhjdhi vanishes. The argument is elementary. Note that dHj/dhj is a function 
of Gjj and Sj. If we now take a derivative with respect to hi when i j, the derivative 
must vanish, because the derivative of Gjj and Sj with respect to hi is zero. Therefore, 



d'^GjiiiUn): „ ,. dT^jdT. 

\h=Q = Gii[lUJn)Gjj[lUJn)Gji[lUJn) 



dhjdhi 



dhj dhi 



\h=0- 



(AlO) 



Evaluating the derivatives explicitly and simplifying the final result then yields Eq. fp7|) . 
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